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IGNITION,  COMBUSTION,  DETONATION,  AND 
QUENCHING  OF  REACTIVE  MIXTURES 


Three  differant  areas  of  the  research  program  have  been  investigated 
during  this  reporting  period 

1  April  1978  through  31  March  1979 
I.  TRANSITION  FROM  DEFLAGRATION  TO  DETONATION 

A.  INTRODUCTION 

The  distance  which  a  deflagration  wave,  traversing  a  combustible  gas 
mixture  in  a  closed  vessel  or  an  unconfined  cloud,  must  travel  before  a 
detonation  wave  is  formed  is  of  great  practical  importance.  Although  in 
some  cases  detonative  combustion  is  desirable  (such  as  in  fuel-air  ex¬ 
plosions,  supersonic  ramjets,  and  high-efficiency  combustion  chambers), 
frequently  it  is  not  wanted  as  in  accidental  fires  of  fuel-air  mixtures, 
and  in  unsteady  combustion  causing  large  pressure  fluctuation  (combustion 
instability)  in  the  combustion  chambers  of  rocket  engines.  For  the  de¬ 
sign  of  devices  which  involve  supersonic  combustion  it  is  extremely  im¬ 
portant  *0  have  criteria  which  can  be  used  to  predict  the  distance  a 
deflagration  wave  must  travel  to  become  a  detonation  wave.  Although  nu¬ 
merous  investigations  have  provided  us  with  a  fairly  good  picture  of  the 
mechanism  of  transition,  there  are  no  expressions  which  permit  us  to  pre¬ 
dict  the  detonation  induction  distance  of  a  combustible  gas  mixture.  The 
research  and  the  experimental  measurements  described  in  the  following 
sections  of  this  report  are  carried  out  to  provide  theoretical  or  empirical 
expressions  for  predicting  induction  distances  in  fuel  oxidizer  mixtures. 

B,  DETONATION,  INDUCTION  DISTANCES  OF  HYDROGEN- OXYGEN- INERT  GAS  MIXTURES 

In  previous  experiments1  with  M.E.R.  mixtures  of  hydrogen  and  oxygen 
(H2  +  \  Oa)  and  hydrogen  +  air  we  observed  that  the  induction  distances 
of  these  mixtures  are  reduced  drastically  when  the  initial  temperature  of 
the  combustible  gas  mixture  is  of  the  order  of  100  K  (see  Table  I).  Con¬ 
sidering  the  fact  chat  the  normal  burning  speed  of  these  mixtures  decreases 
with  decreasing  temperature  it  is  difficult  to  give  the  exact  reason  for 
this  kind  of  a  temperature  effect  on  the  induction  distance.  Examining 
the  various  factors  which  might  affect  the  induction  distance,  we  note 
that  the  density  of  the  gas  mixtures  was  inversely  proportional  to  the 
temperature  in  these  experiments  since  a  pressure  of  1  atm  was  maintained 
in  all  experiments.  Thus  the  thermal  energy  released  per  unit  volume  of 
the  mixture  in  the  low- temperature  mixture  is  greater  than  in  the  mixture 
at  room-temperature ,  Furthermore  the  higher  density  of  the  unburned  gas 
mixture  offers  a  greater  resistance  to  the  expanding  burned  gas.  Because 
of  this  resistance  the  combustion  occurs  more  nearly  at  constant  volume 
than  at  constant  pressure  (no  resistance)  and  consequently  a  rather  marked 
increase  in  pressure  occurs  right  at  the  beginning  of  the  combustion. 

Such  large  and  rapid  pressure  rises  have  indeed  been  observed  at  low 
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Table  I.  Induction  Distances  in  H2  +  £  O2  and 
He  +  +  Qa  +  1.88  Ife  Mixtures  at  300 
and  123  K 


Mixture 

Induction  Distances 

300  K 

123  K 

H2  +  ^  Qa 

H2  +  £  Os  +  1.88  Na 

143  cm 

>600  cm 

15  cm 

230  cm 

initial  temperatur es .  Although  similar  and  even  greater  density  changes 
can  be  produced  by  raising  the  initial  pressure  at  a  fixed  teaperature, 
the  effect  on  the  induction  distance  is  not  as  drastic.  The  results  of 
previous  measurements2  with  hydrogen-oxygen-inert  gas  mixtures, 

H2  +  ^  Qa  +  X,  where  X  ■  Njg,  or  He,  or  A,  are  shown  in  Table  II.  The  ef¬ 
fect  of  pressure  on  the  Ha  +  4  Q2  +  Ife  mixture  is  shown  in  Table  III. 
According  to  these  measurements  the  induction  distances  decrease  only 
moderately  with  increasing  pressure. 

From  theoretical  considerations  it  has  been  found  that  the  induction 
distance  is  a  function  of  the  normal  flame  speed  of  the  combustible  mix¬ 
ture,  the  acceleration  of  the  flame  propagation,  and  the  speed  of  sound 
of  the  unburned  gas.  For  a  given  fuel-oxidizer  mixture  the  flame  front 
accelerates  for  three  reasons;  i.e.,  (1)  the  temperature  of  the  unburned 
gas  is  increased,  (2)  the  motion  of  the  gas  mixture  ahead  of  the  flame 
front  becomes  turbulent,  and  (3)  the  pressure  of  the  unburned  gas  ahead 
of  the  flame  is  increased.  Whereas  the  rate  of  flame  propagation  is  al¬ 
ways  increased  with  increasing  temperature  and  turbulence,  an  increase  in 
pressure  produces  an  increase  in  flame  speed  only  in  fuel-oxygen  mixtures. 
The  flame  speeds  of  fuel -air  mixtures  either  are  practically  independent 
jf  pressure  or  decrease  slightly  at  higher  pressures.  The  role  of  tur¬ 
bulence  is  most  difficult  to  assess.  The  experiments  with  hydrogen-oxygen 
mixtures  at  low  temperatures,  at  which  the  flame  speeds  are  small,  seem 
to  indicate  that  the  role  of  turbulence  is  secondary.  At  the  short  in¬ 
duction  distances  not  much  turbulence  may  have  been  generated. 

Although  f  'une  acceleration  appears  to  have  a  great  effect  on  the 
transition  from  deflagration  to  detonation,  a  more  important  aspect  seems 
to  be  the  rate  at  which  the  pressure  increases  behind  the  flame  front. 

For  instance,  hydrogen-oxygen  mixtures  at  room  temperature  have  rather 
large  induction  distances,  whereas  acetylene-oxygen  mixtures  whose  burning 
velocities  are  practically  the  same  have  much  shorter  induction  distances. 
The  major  difference  between  the  two  mixtures  is  the  density.  .The  density 
of  an  M.E.R.  acetylene-oxygen  mixture  is  2.5  times  as  great  as  that  of  an 
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Table  XI.  Induction  Distances  in  Ha  +  £  Oa  +  nX  Mix 
tures  at  1  and  5  atm  Initial  Pressure  and 
Ti  *  313. J-5  K  (Length  of  Tube  -  350  cm; 
I.D,  *  5 


nX 

Pi 

?ind 

(atm) 

(cmj 

£  NS 

1 

5 

230 

120 

1  Ife 

l 

5 

>350 

250 

1  He 

1 

180 

5 

70 

£  He 

1 

5 

0  0 

£A 

1 

5 

180 

45 

1  A 

1 

5 

240 

115 

£  CQs 

l 

5 

240 

115 

Table  III.  Induction  Distance  as  a  Function  of  Pressure  (Ti 


Pressure 


0.5  atm 


Mixture 


H2  +  £  02  +  N; 


5  cm 


325  cm 


270  cm 


M.E.R,  hydro  gen-oxygen  mixture.  Although  the  combustion  enthalpy  of  a 
unit  volume  of  the  acetylene-oxygen  mixture  is  also  much  greater  than 
that  of  the  hydrogen- oxygen  mixture,  the  difference  may  not  be  significant 
since  both  mixtures  have  practically  the  same  flame  temperature  so  that 
a  constant  volume  combustion  produces  very  nearly  the  same  pressure  ratio 
of  <-10  when  the  initial  gas  temperature  is  300  K.  The  formation  01  :dgher 
pressures  during  the  initiation  of  the  detonation  wave  results  from  an 
energy  transfer  because  work  is  done  by  the  burnt  gas  on  the  gas  behind 
the  shock  wave  preceding  the  flame.  The  role  of  density  is  also  evident 
when  we  compare  the  induction  distances  of  H2  +  £  Q2  +  nA  mixtures  with 
those  of  He  +  ■£■  O2  +  nHe  mixtures,  as  shown  in  Table  II.  The  energy 
density  in  all  three  mixtures  is  the  same.  The  burning  speeds  of 
Ha  =  £  O2  +  nHe  are  much  greater  than  those  in  Ha  +  ^  Qa  +  nA  mixtures, 
whereas  those  in  the  mixture  containing  nitrogen  are  slightly  smaller 
(Fig.  1). 

The  induction  distances  with  argon  are  much  smaller  than  those  in 
the  mixtures  with  nitrogen  or  helium.  Diat  the  induction  distance  in 
mixtures  containing  nitrogen  is  practically  identical  with  that  in  the 
mixture  with  helium  can  be  explained  on  the  basis  of  the  lower  adiabatic 
combustion  as  well  as  detonation  temperature  of  the  mixture  containing 
nitrogen.  The  large  induction  distance  of  Ha  +  ^  Qg  +  CQ2  mixtures  in 
spite  of  their  high  density  must  be  attributed  to  the  interaction  between 
COa,  Ha,  and  H2O.  For  the  same  reason  the  flame  speed  and  flame  tempera¬ 
ture  of  this  mixture  are  also  very  low.  Experiments  with  these  mixtures 
at  low-temperature  will  be  carried  out  in  the  near  future. 

C.  FLAME  SPEEDS  OF  HYDROGEN-OXYGEN- INERT  GAS  MIXTURES 

To  determine  the  qualitative  and  quantitative  relationships  be Ween 
the  flame  parameters  and  the  induction  distances  of  various  combustible 
gas  mixtures,  flame  speeds  are  being  measured  for  several  hydrogen- oxygen- 
inert  gas  mixtures  with  different  amounts  of  additive  gas  and  over  ,s 
wide  a  temperature  range  as  possible  (100  to  400  K  provided  condense. tion 
of  the  additive  gas  does  not  occur) .  A  special  nozzle  burner  wit.  a  sine 
curve  contour  is  used  for  these  measurements  to  obtain  laminar  flames 
with  nearly  straight  flame  cones  even  at  Reynolds  numbers  up  to  10,000. 

The  measurements  are  extremely  reproducible.  Results  obtained  so  far  are 
shown  in  Fig.  1.  It  is  interesting  to  note  that  H2  +  -|  Q2  +  argon  mix¬ 
tures  have  practically  the  same  flame  speeds  as  H2  +  ■§•  Q2  +  nitrogen  mix¬ 
tures  in  spite  of  the  large  physical  and  chemical  differences  between 
these  two  additives. 

D.  EFFECT  OF  INITIAL  TEMPERATURE  ON  WAVE  SPEED  OF,  AND  FINAL  PROPERTIES 
BEHIND  C.J.  DETONATION  WAVES  IN  THE  M.E.R.  MIXTURE  OF  HYDROGEN  WITH 
AIR 

The  wave  speed,  gas  speed,  pressure,  temperature,  and  density  were 
calculated  for  the  M.E.R.  (maximum  energy  release;  molar  ratio  of  hydro¬ 
gen  to  oxygen  =  2:1)  mixture  of  hydrogen  and  air  at  various  initial 


temperatures  (Ti  w  100  to  473  K)  and  initial  pressures  (pj.  *  1  to  100 
atm).,  The  results  cure  shown  in  Pigs.  2,  3,  4,  5>  and  6.  According  to 
these  calculations  the  wave  speed  !:•  inverse.!;''  proportional  to  the  tan- 
perature  of  the  unburned  gas  but  proportional  to  its  pressure  (Fig.  2), 
whereas  the  gas  speed  increases  very  slightly  with  increasing  initial 
temperature  and  pressure  (Fig.  3)  .  Although  the  final  temperature  in¬ 
creases  with  increasing  initial  temperature,  the  increase  is  much  less 
than  that  of  the  initial  gas  temperature.  This  behavior  is  probably  due 
to  the  less  efficient  energy  transfer  from  burned  gas  to  the  gas  behind 
the  wave  which  will  be  discussed  in  the  next  section.  The  increase  in 
Tj  caused  by  an  increase  in  pL  is  due  to  a  reduction  in  dissociation  of 
the  combustion  gas  molecules  (Fig.  4) .  Because  of  the  importance  of  tur¬ 
bulence  on  the  rate  of  flame  propagation  the  effect  of  pressure  on  the 
speed  of  the  detonated  gas  is  shown  in  Fig,  5,  which  clearly  reveals  that 
u3  increases  only  very  slightly  with  pressure  although  the  induction  dis¬ 
tances  of  almost  all  combusti  ble  mixtures  decrease  rather  markedly  with 
a  rise  in  pressure  of  the  initial  gas.  The  most  significant  observation 
is  the  large  increase  of  the  pressure  of  the  detonated  gas  as  the  tem¬ 
perature  is  lowered  (Fig.  6).  On  the  other  hand,  a  change  of  the  initial 
pressure  hardly  affects  the  pressure  ratio.  As  shown  in  the  next  section 
the  relative  amount  of  energy  transferred  from  the  downstream  burned  gas 
to  the  gas  behind  the  normal  shock  wave  preceding  the  combustion  zone  is 
directly  proportional  to  the  pressure  ratio;  e.g., 


where  Rj.  =  . 

Although  it  is  tempting  to  relate  the  drastic  reductions  of  the  induction 
distances  at  very  low  initial  temperature  of  the  unburned  gas  to  the  en¬ 
ergy  transfer,  no  definite  statement  can  be  made  at  this  time  because  the 
observed  pressures  of  the  detonated  gas  mixture  at  low  initial  tempera¬ 
tures  were  found  to  be  much  lower  than  the  calculated  values.  Therefore, 
more  experiments  with  different  combustible  gas  mixtures  at  low  initial 
temperatures  are  necessary  to  establish  a  quantitative  relationship  be¬ 
tween  mixture  properties  and  detonation  induction  distance. 


XI.  INVESTIGATION  OF  THE  MECHANISM  OF  FLAME  QUENCHING 
The  following  topics  were  investigated: 

(1)  The  effect  of  salt  coatings- -thick  coatings  of  salt  were  obtained  by 
painting  solutions  of  potassium  chloride,  sodium  bicarbonate,  and 
tribase  potassium  phosphate  on  both  surfaces  of  the  quenching  plates. 
Thick  layers  (about  0.16  cm)  of  the  salts  were  used  so  that  the 
flame  could  react  chemically  with  the  salt  during  the  prolonged  con¬ 
tact  of  the  flame  with  the  quenching  plate  surfaces.  It  was  observed 
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that  salt  coatings  of  potassium  chloride  and  sodium  bicarbonate 
vanish  rather  rapidly  when  exposed  to  the  flame,  whereas  the  coating 
of  trlbase  potassium  phosphate  remains  on  the  plate  surface  much 
longer.  It  was  found  that  the  appearance  of  the  flames  of  methane- 
air,  methane-oxygen,  acetylene -air,  and  hydrogen-air  mixtures  was 
affected  by  the  three  salts  although  the  quenching  distances  were 
hardly  changed.  Therefore,  it  can  be  concluded  that  these  salt 
coatings  do  not  promote  the  flame  quenching  process. 

(2)  The  effect  of  gas  speed  on  quenching  distance — further  experiments 
were  performed,  by  increasing  the  volume  flow  rate  of  the  unburned 
gas  to  200  cc/s  (twice  the  volume  flow  rate  used  in  previous  experi¬ 
ments).  For  the  gas  mixtures  considered  (i.e.,  methane-air,  methane- 
oxygen,  acetylene-air,  and  hydrogen- air ) ,  the  quenching  distances 
remain  unchanged.  Consequently,  it  can  be  concluded  that  for  laminar 
flames  the  measured  quenching  distances  are  independent  of  the  ini¬ 
tial  gas  speed. 

(3)  The  mathematical  expression  developed  for  the  flame  quenching  dis¬ 
tance,  dq;  i.e.. 


4Xh  +  2Dmlxcpb 


W 


where 


sugj 


*  thermal  conductivity  of  the  gas  mixture  at  average  flame 
temperature, 

■  height  of  the  quenching  plate, 

=  diffusion  coefficient  of  the  gas  mixture  at  average  flame 
temperature, 

«  specific  heat  at  constant  pressure  at  an  average  flame  tem¬ 
perature, 

«  short  side  of  the  burner  port, 

=  flame  speed,  and 

*  geometrical  factor  of  the  order  of  unity, 

sts  that  the  quenching  distance  should  be  affected  by 

(a)  the  width  b,  of  the  burner  port  and 

(b)  the  height  of  the  quenching  plate,  h. 


In  order  to  study  the  effect  of  the  width  of  the  burner  opening, 
experiments  were  performed  with  methane -air  flames  by  using  a  new 
burner  whose  width  was  two  tines  that  of  the  original  burner;  i.e., 
0.3*4  cm.  the  measurements  reveal  that  the  quenching  distance  is 
slightly  reduced  (of  the  order  of  10  percent)  when  b  is  increased. 

Then  experiments  were  performed  with  the  original  burners  after  the 
height  of  the  quenching  plates,  h,  was  increased  two  times  the  orig¬ 
inal  height  (from  2.5  to  5  cm).  These  results  showed  that  for  the 
gas  mixtures  examined,  the  quenching  distances  are  independent  of 
the  quenching  plate  height,  h. 

(4)  In  order  to  determine  the  relationship  between  the  quenching  dis¬ 
tance  measured  in  our  laboratory  and  those  established  by  an  entirely 
different  technique  by  previous  investigators,  experiments  were  per¬ 
formed  to  measure  that  distance  between  the  quenching  plates  (once 
the  quenching  is  established)  where  the  flame  burning  on  the  top  of 
these  plates  just  settles  down  suddenly;  i.e.,  flashes  back.  To 
determine  this  point  the  quenching  plates  (made  of  copper)  were  moved 
apart  very  slowly.  Sudden  flash-back,  however,  did  not  occur;  the 
flames  entered  gradually  and  slowly  descended  between  the  two  plates. 
Further  studies  are  necessary,  probably  with  temperature-controlled 
plates,  to  explain  this  unexpected  behavior. 

(5)  Further  experiments  were  performed  by  closing  the  side  openings  of 
the  copper  plates.  In  this  way  a  rectangular  enclosure  of  the 
quenching  zone  was  formed.  It  was  found  that  at  a  certain  di stance 
between  the  plates  the  flames  burning  at  the  top  of  these  plates 
suddenly  jumped  down  to  the  burner  tip  when  the  plates  were  moved 
apart  slowly.  When  copper  plates  were  used  all  around,  it  was  not 
possible  to  see  where  the  flame  was  quenched  and  furthermore  it  was 
not  possible  to  see  whether  the  quenching  process  was  sudden.  In 
order  to  observe  the  flash-back  point,  glass  plates  were  employed  to 
form  the  enclosure.  It  was  found  that  the  flame  entered  the  burner 
suddenly  and  when  the  plates  were  moved  together,  at  a  certain  dis¬ 
tance  the  flame  jumped  to  the  top  of  the  plates.  The  distance  be¬ 
tween  the  plates  at  which  the  flame  suddenly  leaves  the  burner  and 
jumps  to  the  top  of  the  plates  was  measured  for  methane-air  flames, 
the  results  are  compared  below. 
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i 

% 

Methane 

Quenching  Distance 
without  Glass  Plates 
(am) 

Quenching  Distance 
with  Glass  Plate 
Enclosure 
(am) 

7.0 

2.43 

3.68 

8.0 

2.22 

3.28 

9.5 

2.10 

3.16 

11.0 

2.25 

3.72 

12.0 

2.41 

5.30 

It  mast  be  concluded  that  the  true  quenching  distance  Is  given  by 
the  values  obtained  with  the  enclosure  of  the  narrow  side  of  the 
quenching  space. 


in.  NEW  CALCULATION  TECHNIQUES 

A.  CALCULATION  OF  DETONATION  PARAMETERS  AND  ENERGY  TRANSFER 

For  the  present  calculations  we  consider  the  mixtures 

Ha  +  £  Qa  +  v  Y 
¥ 

where  Y  *  argon,  helium,  nitrogen,  or  carbon- dioxide.  For  estimated  val¬ 
ues  of  T3  and  p3  and  specified  values  of  Ti,  Pi,  v,,  and  Y  the  qt>3 

^HaO*  ^Qa*  ^Ha*  V  V  ^He*  ^A*  ^Na*  T,N0»  ^CQa*  ^  W  are  cal” 

culated  according  to  the  technique  explained  In  Section  III-B. 

Since  the  adiabatic  flame  temperature  at  constant  volume  (t3v*^ )  is 
the  lowest  temperature  a  supersonic  combustion  wave  can  attain,  it  is  ad¬ 
visable  to  calculate  this  temperature  first.  In  this  case  we  have 
wx  *  w3  «  «  and  u3  »  0.  To  find  the  pressure,  p3,  which  is  compatible 
with  the  assumed  tenperature,  T®  ,  we  use  the  condition  that 

pgalc  .  Pl  .  »  when  ■  1  , 

where  1fl3  «  ]T  Tjt>3Wt  . 


Since  for  a  given  T3  the  molecular  mass  7>13  increases  as  p3  is  increased 
(less  dissociation)  and  vice  versa,  pS^0  will  be  aaaller  than  the  correct 


value  when  pSst  is  larger  than  the  correct  value  and  vice  versa.  Hence 
a  subsequent  estimate  of  p3  is  obtained  from  the  interpolation  formula 

est(n)  n  est  .calc 

P3  *  (i  -  s;  P3  +  6P3  » 

where  0  <  ft  <  1  ,  usually  0.5. 

The  iterations  are  repeated  until 

calc(n)  est(n) 

100  •  ra  "juaS —  <  ■  o-001*  • 

P3 

After  the  correct  value  of  p3  has  been  calculated,  the  energy  equation  is 
used  to  calculate  the  correct  temperature  T3.  For  combustion  at  constant 
volume  ve  have 


which  can  be  written  in  terms  of  the  reduced  formation  enthalpies  as 
follows:  first  we  write 

^2  tef .  5.  tef  , 

^  \RT/  W  \OT/ 

•solving'  for  T3  we  obtain 


Since  usually  only  formation  enthalpies  are  tabulated  we  use  the  relation- 
ships  e  «  h  -  pv  and  E  ■  H  -  3CT  (thermally  perfect  gas)  and  obtain 
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and 


/HA1  /H-Eo\ 

Vsf/,  •  w, 


AEq,  i 


The  ^Eo»i  *  Eo>  i  -  £  v^*'Eo,el  may  considered  as  the  formation  ener- 

el  calc  est 

gies  (or  enthalpies)  of  the  species  i  at  T  »  OK.  When  T3  j  13  the 

calculations  are  repeated  with  improved  estimates  according  to  the  equa¬ 
tion 


Test(n) 


T|St(n-l)  _  L 


where 


until 


TestO 


100  . 


r  _est(n-l)  1 

/hf\ 

Ti  /: 

hf\Tl 

l\W  '  \ 

■  wlv 

‘  \ 

Test(n)  _  Tcalc(n) 


Test(n) 


<  0.001$  . 


Now  we  can  proceed  to  calculate  detonative  combustion  parameters. 

For  assumed  values  ot  T3(T3  >  Tfu*°)  and  p3(p3  >  ps*,*°)  first  we  must  cal¬ 
culate  again  the  rji>3.  To  find  the  correct  pressure  for  the  assumed  tem¬ 
perature  we  'solve'  for  p3  from  the  energy  equation  which  yields 


_calc 

Pa 


T3  trij.  /hf\T3  fa 

Ti  *  W3  \RT  /  ~  \RT 

1  +  Sa  .  El  .  J!l 

Ti  P3  TMa 


hf\Tl 


If  pi^0  ft  p3st,  the  calculations  are  repeated  with  improved  estimates. 


est(n)  ,x  est(n-l)  .  calc(n-l) 

Pa  v  =  (1  -  5)  Pa  v  '  +  6P3  , 

I  est(n)  calc(n) I 


until 


caic(n) 

P3 


<  £sp3  =  0.0001$  . 


Any  compatible  pair  of  T3  and  p3  values  represents  a  solution  of  the 
Hugoniot  equation  and  thus  a  supersonic  combustion  wave.  We  have  a  weak 


detonation  when  the  selected  temperature,  T3,  is  greater  than  T3  hut 

C  J 

smaller  or  just  equal  to  T3  ,  the  temperature  of  the  so-called  Chapman  - 
Jougnet  detonation  wave  which  is  characterized  by  the  condition  *3  =  wa>3. 

C  J 

Weak  detonation  waves  involve  no  shock  transitions.  When  T3  >  T3  we 
have  a  strong  detonation  wave  in  which  a  subsonic  combustion  wave  occurs 
behind  a  normal  shock  wave  which  traverses  the  unburned  gas  mixture  at 
the  speed  w.  For  strong  detonation  waves  we  have  W3  <  wa,3. 

To  calculate  the  wave  speed  and  the  conditions  behind  the  wave  of  a 
Chapman- Jougnet  detonation,  which  is  the  one  which  occurs  as  a  traveling 
wave,  we  meke  use  of  the  condition  that  at  the  tail  of  the  wave  the  wave 
speed,  W3,  is  equal  to  the  speed  of  sound  in  this  gas,  wa,3.  This  condi¬ 
tion  can  be  used  to  develop  the  following  equation  for  a  calculated  value 
of  T3: 


b/.J  -  Av 


These  parameters  must  be  calculated  for  the  estimated  temperature,  T3  . 

C&lC  i  est 

When  T3  f  T3  the  calculations  are  repeated  with  an  improved  estimate 
which  is  obtained  from  the  following  expression: 


T|«(n) 


Test(n-1) 


where 


Test(n-1) 


est(n-l)  m 

\RT/  W  *  \RT/  * 


The  iterations  are  continued  until 


^est(n)  _  Tcalc(n) 


<  0.001^  . 


After  the  correct  temperature  has  been  determined  the  speed  of  the  deto¬ 
nation  wave  is  calculated  by  means  of  the  following  equation: 


T  Ea  _  1  f 

!.J.  .  gL  .  PI  M. 

w  nrnrrs: 

Ti  p3  yn  3_ 


Although  these  calculations  produce  data  which  describe  only  the 
steady-state  motion  of  the  wave  and,  therefore,  cannot  be  used  directly 
to  calculate  the  induction  distance,  it  appears  that  the  pressure  ratio, 

and  the  energy  transfer  from  the  burned  gas  to  the  gas  behind  the  wave 

may  be  related  to  the  magnitude  of  the  induction  distance  as  determined 
in  these  calculations.  The  amount  of  energy  transfer  can  be  calculated 
as  follows:  The  stagnation  enthalpy  of  the  unburned  gas  relative  to  a 

laboratory-fixed  coordinate  system  is  h^1  and  that  of  the  gas  at  the  tail 

of  the  wave  is  hTa  +  Hence  the  amount  of  energy  transferred  is 


Ah'"  -  hTa  +  ^  -  hTl  =  hf3  + 

Relative  to  the  wave  we  have 

v,T3  wf  .  Ti  w£ 
hf>  +  2  *  hf  +  2 


-  hf1  . 


and  thus 
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Substitution  leads  to 


.  w£_|_sS  .  ±*A.  „lW3  . 


W1W3 


or 


'  V1  '• 


so  that 


^;  =  E2. 1 
m\  pi  x 


This  expression  shows  that  the. relative  Energy  transfer. 


£h 


,  is 


proportional  to  the  pressure  ratio,  .  The  relationship  between  these 

parameters  and  the  induction  distance  will  be  examined  in  future  experi¬ 
ments. 


B.  CALCULATION  OF  PERFORMANCE  OF  SUPERSONIC  RAMJET  ENGINE 

The  analysis  includes  the  following  processes  and  engine  parameters: 

1.  Diffusion  from  the  Free-Streaii  Air  Speed,  u  ,  to  a  Specified  Value, 
Upjj,,  at  the  Diffuser  Exit  or  the  Combustion” Chamber  Inlet 

(Uc,l  =  ude)* 

a.  The  diffuser  is  assumed  to  be  isentropic. 

1.  The  decelerated  airflow  is  subsonic 

2.  The  decelerated  airflow  is  supersonic 

b.  A  normal  shock  wave  is  formed  at  the  inlet  of  the  diffuser  and 
subsonic  diffusion  to  u^  follows. 

c.  Hie  deceleration  of  the  incoming  airflow  occurs  first  via  an 
oblique  shock  wave  on  a  wedge,  then  via  a  normal  shock  wave  which 
is  followed  by  a  subsonic  diffusion  to  u^.. 
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2.  Fuel  Is  Added  to  the  Air  Stream  Tangentially  and  Burned  In  a  Constant 

Are*  Ao  *  aBE  +  *1 UEL' 

a.  The  airflow  entering  the  combustion  chamber  Is  subsonic. 

b.  The  airflow  entering  the  combustion  chamber  is  supersonic, 

(1)  Combustion  occurs  without  any  shock  waves  ( stationary  weak 
detonation  wave) 

(2)  Combustion  is  preceded  by  a  normal  shock  wave  (stationary 
strong  detonation  wave) 

3.  Expansion  Qhrougi  Exhaust  Nozzle 

a.  The  expansion  is  isentropic:  convergent -divergent  nozzle  with 
Se  "  *a<Me  >  15  • 

b.  Die  nozzle  consists  of  a  convergent  section  only  (M  '  =  1)  p  >  p  . 

6  6  a 

Flow-through  nozzle  is  assumed  to  be  isentropic  but  expansion 
from  p  top  outside  the  nozzle  is  accompanied  by  a  rather  large 

increase  in  entropy. 

4.  Calculation  of  Important  Cross-Sectional  Areas ;  Only  Relative  Values, 
Referenced  to  the  Inlet  Area  Ai,  Will  be  Given. 

a.  Apj.  =  throat  area  of  diffuser  for  the  case  of  isentropic  diffusion 
to  a  subsonic  speed  u_p. 
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tn.,_  *  28. 85297479  kg/kmol 


air 
FUH,:  Ofe. 


Ha 


v„  -  (to  be  specified 

Jk 


m 


I-'OEL  *  ^ 12,01  +  vHa  *  2*01^)  4g/lcnol  s  tTTCHg  ^ 


FUEL-AiR  MIXTURE:  CH^  +  *  3.76  Na) 

He  '  1  ■ 


Ha 


MIXTURE  RATIO:  f 


■WVr  -  WeB2vJ! 


^  vQa  -*W 


Vq^  *  to  be  specified 


All  underlined  properties  and  At  ^or  must  be  specified. 


Furthermore,  the  thermodynamic  properties  of  the  fuel,  the  air 
and  all  combustion  products  must  be  available. 

Tt  «  integer  hundred  nearest  T  but  T-  <  T. 

Li  U 

1  (1't  - T.)/100 
T 

c  .  ■  Entropy  Coefficient  of  Species  i  at  Temperature  T. 
s,i 

For  air  at  T  we  have 
00 

V 
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iTJ 


The  depend  on  and  p^.  They  must  be  calculated  by  an 
iterative  procedure.  For  estimated  values  of  and  ^ob¬ 
tained  from  Eq.  (2)  and  Eq.  (3)  with  estimates  of 
(h^RT^gE^DE  and  c^DE;  approximately  hf/RT  ~  7/7  -  1  and 
c  ~  7/7  -  1)  the  r] .  are  calculated  as  follows: 

S  /  1  y  Dili 


STEP:  (1) 


est 

—COW  _ 

aDE  L 


(2)  [a(0)(TL)](1  “DE)  .  [a<°)(IL  ♦  100)]"1*  .  . 


•  (O'  - *(0) 

(3)  .  r«<H0>/T.  ♦  lOO^"1®  . 


K 


<7>  "so  *  •  Jc  - 

1(fe  *  (life  +  5  1no)  *  3-76  "  (\  +  I'no)  *  2  • 

When  tj^10  t  repeat  Steps  (5)  through  (8)  with 

(9)  i“t(n)  -  s  •  n£te<»-«  +  0--5)  -  , 

where  0  <  5  <  1  (usually  l  »  0.5  but  in  some  cases  5  <  0.1  to 
obtain  convergence).  Repeat  iteration  until 

^calcCn)  _  ^est(n) 

<10>  100  •  ■  cali'fn?8 -  <  ‘  O-00001*  • 

^oa 

Then  calculate 


(H)  WDE  =  •  (2^  +  n0  +  t)nq)  ►  Wjj  •  (2rjNa  +  t]noj 
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“  PDE  *  PDE  »  repeat  Steps  (2)  through  (14)  with 

(15)  p*f tn+1)  -  S  •  P°fc(n)  *  (1  -  !)  •  P^t(n)  , 

where  i  =  0.5.  Continue  the  iteration  until 


(16)  100  . 


calc(n)  est(n) 


calc(n) 


<  £p_p  =  0.004  . 


Then  calculate 
/hf\TDE 

<i7>  U)  ■ 


(mL  •  (i  -  -w  *  m) 


Tr+100 

1i 


>o  4  (te8)^  •  (x  -  "Be)  *  (^)c 

*  •  to!1  •  (x  -  "be)  ♦  m 

L  "Z 


TL+100 


Tl+100 


■  (>  -  w  *  (*;■”  -vH 


When  +  T*^  ,  repeat  steps  (1)  (see  p.  23)  through  (18)  with 

(19)  -  8T“1C<n>  Ml  -  I)  •  !«*<“>  , 

where  |  -  0.5.  Continue  the  iterations  until 
Tcalc(n)  Test(n) 

<*»  100  •  ~  SSSX& <  "be  -  • 

ade 


The  basic  configurations  of  the  engine  may  be  pictured  as  follows: 


l.b.  Normal  Shock  Wave  at  Inlet 


Instead  of  Eqs.  (2)  and  (3)  (see  p.  22)  use  now 


and 


EST 

p2 


VC  -  (y  '  1) 

7  +  1 


let  7  ~  1.3 

for  M  ~  10 


(4) 


(5) 


for  the  estimates  of  temperature  and  pressure  behind  the  normal 
shock  wave  attached  to  the  inlet.  Then  use  steps  (1)  (see  p.  23) 
through  (11)  to  calculate  the  t]t>2  (steps  1-11  are  cannon  to  both 
calculation  procedures) .  Then  calculate 


(12*) 


n 


NO 


(1  -  *2) 
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This  calculation  of  gj-  is  Identical  to  that  of  8t«p  (17)  of  the 

previous  procedure  with  the  exception  that  Tnip  has  been  replaced 
by  Ta.  Then  calculate  ^ 


(13*)  P l*C 


-  P.  *  B  +  ^B®  + 


3F 


tt! 


2  '"air 

TT  *  “W“ 


where 


B  * 


**“  Walt/M1?1  /“ff.  i/*f*  V 


Gft 


-1)- 


When  pa*10  t  p|St»  repeat  steps  (1)  through  (11),  (12*),  and  (13*) 
with 

(14*)  p|st^a+1)  -  {pi*10^  +  (1  -  0  pf8tW  , 
where  5  »  0.5.  Continue  the  iteration  until 


(15*)  100 


calc(n)  est(n) 

*  -  nSf — 


calc( 
Pa 


<  AP2  =  0.001* 


Then  calculate 


(16*)  Tg*10  »  T 


/>,  \T°° 

UfJ  +  °'7h 

XRT/air 


iy£  U  - 


„EST 


m 


’air 


TO**  1 

Wa  P2 


air 


/MTf1 

WM 


When  ^  T^fst,  repeat  steps  (1)  through  (11)  and  (12*) 

through  (16*) 

(i7»)  iS5**"*11  .  s  .  il*1*  +  (i .  o  i|5t , 

where  5  is  as  before.  Continue  the  iteration  until 


(18*)  100 


Tcalc(n)  _  Test(n) 
~S  T|alc(n) 


<  <^r2  =  o.ooi*  . 
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Then  calculate 


x  $*)  u2  "  u« 


Ta  »-  n* 

*1  '  *2  ‘  ~w  ■ 


If  ua  la  larger  than  practical,  subsonic  diffusion  from  ua  to 

is  necessary,  These  calculations  are  identical  with  those  given 
in  Section  l.a.  when  the  subscript  «  is  replaced  by  the  sub¬ 
script  2. 

l.o.  One  Oblique  Shock  Wave  Followed  by  a  Normal  Shock  Wave 


Use  the  equations  of  Section  l.b.  but  replace  u^  by 

u  ,  »  u  •  sin  a,  where  cr  is  the  angle  of  the  shock  wave  formed 

by  the  wedge  (see  sketch  above;  cr  is  given ! ) .  Replace  1^  by 

M  _  »  M  •  sin  <r.  With  the  values  of  u  (M  ) ,  T  ,  and  p 
°>,n  »  «,n  «,n'  7  a  -“a 

calculate  the  state  of  the  air  behind  the  oblique  shock  by  the 

procedure  outlined  in  Section  l.b.  These  calculations  lead  to  pa, 

Tg,  and  u0  .  The  gas  speed  parallel  to  the  wedge  is 
n 

Ug  *  jr;  +  (Uw  .  cos  <r)2  (Mg  >  1) 
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and  the  half  angle  of  the  wedge  Is 

&  =  <7  -  tan"1  . 

\  « 


COS  <J\ 


Then  calculate  the  conditions  behind  the  normal  shock  wave  which 
are  designated  by  the  subscript  3.  For  the  first  estimate  of  T3 
and  P3  use  instead  of  Eqs.  (4  and  5): 


*P*  -  b 


(z~P)  •  a  •  ^ 


EST 

P3 


- 


Iuse 

7  -  1.: 


If  Ug  >  UpE,  it  is  necessary  to  employ  subsonic  diffusion  from 
u^  to  uDE  by  means  of  the  calculations  developed  in  Section 
l.a.  The  subscript  «  is  replaced  by  the  subscript  3.  In 
step  (16*)  the  term  0.7MI,  is  replaced  by  *  W2. 


2.  Calculation  of  State  of  Combustion  Gas 
From  continuity  and  momentum  equations 


[■»<•«>  -  u  1  , 
L  c,max  cj  * 


where 


(M-0)  . 

c,max 


*  *  UFUEL 

,  .  PFJEL 

PFUELUflJEL. 

+  UDE 

r  +  it 

L1  +  ’™ 

l  1  +  f  ) 

From  the  energy  equation  we  obtain 


\  c,max/  2  TFT  \RT/n  ’ 


where 


Substitution  of  Eq.  (8)  into  Eq.  (9)  leads  to 


Since  tj.  depend  on  T  and  p  ,  the  calculation  of  the  state  of  the 

c  c 

combustion  gas  u^,  consists  of  an  iterative  procedure  comprising 

three  loops.  The  inner  loop  consists  of  the  calculation  of  the  . 

Use  T?ST  ~  T«_  +  1500  K  and 
c  DE 
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w/ 


yW 

C8  *  (^HaO  +  ’’Ha  +  I  ’’OH  +  I  ’’h)  *  ’’CQa 
Vife 

"  2  (’’CO  +  ’’feO  +  ’’OH  +  ’’o  +  ’’no)  * 
when  T]?®10  ft  repeat  calculations 

M2  M2 

with 

n“T<ntl)  -  6^c(n)  *  (i  -  5) 


In  the  second  loop  the  p  is  calculated  so  that  it  is  compatible  with 

the  estimated  temperature  From  the  continuity  equation  and  with  the 

condition  that  c 


5k 


Ac  “  +  AFUEL 


we  obtain 


¥  TT  ’ 

c  c 


_ 1 4-  f 

t  TPE 

apFUELUFUEL  ^1DEUDEPDE 


When  p®8,10  ^  p®3^  all  previous  calculations  are  repeated  with 

pest(n)  mi  .  pcalc(n-l)  +  (1  _  ft)  pest(n-l) 


until 


calc(n)  _  est(n) 
*c  *c 

_calc(n) 


<  &p  =  0.0001$  . 

c 


Now  we  calculate  the  formation  enthalpy  of  the  combustion  gas; 


'h-\TC 


"02  -  ».)  ♦ 

(<*£  (>•■.)•  w; 

V*  ((tI6)^  (l  •  "c)  '  (t^) 


Tt +100 

,  ii 


Tt  4-100 

i 


Tt  4-100 
i* 


•  m  1 4- 

7 


•  m  i  4- 


•  m  1  4- 


,V100  47286 1 

CCfe  '  c  '  Tf1, 

c 


(»£<*  --.)*»£ 
>  (i  ■  ■•)  •  (») 


T_  4-100  T  3  Aft  A  \ 

CO  c  jEST  I 


TL+10°  287^6 \ 

HsO 
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'oh  i1  ■  m=)  +  •  “o  +  !|S§)  + 


•  mc  +  mEST 


and  then  u  according  to  Eq.  (10)  and  T  according  to  Eq.  (8) .  When 
c  c 

Tc&lc  ft  Test  all  iterations  (r\.  and  p  \  have  to  be  repeated  with 

C  C  \  1,C  C/ 

■jESTCn)  ^  gTcalc(n-l)  +  (X  .  ftjest(n-l) 

.  calc  ^  1  (M-C) 

when  uc  <  g  , 


^EST(n)  =  (1  +  d  TEST(n-l)  _  |Tcalc(n-l)  ^ 


wh«n  >  |  u<»*> 

c  2  c,max 

3.  Calculation  of  Exhaust  Speed,  u  ,  Specific  Thrust,  and  Thrust 
Specific  Fuel  Consumption  e 

a.  The  Expansion  is  Assumed  to  be  Isentropic  and  Complete  So  That 
Pe  -  pa 


With  an  estimate  of  the  temperature  at  the  nozzle  exit  Tram  the  sim¬ 
plified  entropy  equation  the  q,  are  calculated  by  the  same  scheme 

i,e 

as  used  for  the  q,  calculation.  Then  the  nozzle  exit  temperature 

i,c 

is  calculated  from  the  entropy  equation 
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'4* 


i°Uc  ■  *i8» 

e  c 


where 


,X. 


£  1i,< 


Z'Kk  ».  i 


LPc 


‘Of 


e  L  *«  T  J 


Co-LC  /  AST 

When  T  #  T  ,  the  calculations  are  repeated  with  new  estimates  of 
e  e 

Tg  obtained  from  the  equation 

Test(n)  x  calc(n-l)  Test(n-1) 

e  5  e  '  s/  e 


until 


100 


_calc(n)  _  Test(n) 

e _  e 

mCalc(n) 


<  0.001$  . 


The  speed  of  the  exhaust  jet  is  then 


ue  "  \R  *  a'e 


-  2hJe 

f,c  f,e 


or 


(l 

if  \t=  «.  ii 

if\Te 

tr/c  ~^v? 

®ij 

ue  =  \/uc  +  2 


The  specific  thrust  is 

Fa  '  (1  +  f)  ue  -  “f 

and  the  thrust  specific  fuel  consumption  is 


sfc 


=  3600  i  r^i 


L  N  J  * 
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b.  The  Expansion  is  Incomplete  (Convergent  Nozzle  with  M  ■  1) 

For  an  estimated  value  of  the  exit  temperature  ^Tft  ■  ■  0.87  Tfi^ 

together  with  an  estimate  of  the  pressure  at  the  nozzle  exit, 

are  calculated  first.  Then  the  entropy 
the  pressure  which  is  compatible  with 

_ 

•fir  E  -  E  it,.*")!,. 

e  c 

rn  '  w 

Ac 

When  pca^c  4  pest  the  iteration  is  continued  with 
0  0 

pest(n)  =  5pcalc(n-l)  +  ^  ^  pest(n-l) 

0Q  0 


e 


p' 


JSBX 


0.53  Pc, 


the  tj 


i*e 


equation  is  used  to  determine 
the  estimated  temperature 


calc 


K 


Icse 

e 


er 


A 


until 


100 


pcalc(n)  _  pest(n) 


P, 


calc(n) 


<  4>e  «  0.001#  . 


^Note :  for  each  new  p®st  a  new  set  of  g  must  be  calculated.^ 
Now  ,h  the  energy  equation 

2  2  T  ^e 

and  the  condition  that  uf  =  wf  =  7Ze  x=r  > 

6  8»)  6  0  me 


a  calculate!? 

f ..’..‘Si  of  the 


T  = 
e 


it  temperature,  T  ,  is  determined  from  the  following 
vgy  equation: 


which  leads  to 


When  T^8^0  /  T88^,  all  three  Iterations  (r\.  ,  p  ,  and  T  \  are  re 

0  ©  \  0  0  / 

peated  with 

Test(n)  =  iTeale(n-l)  +  ^  ^  Test(n-1) 

0  0  0 


Then  the  exhaust  speed  is  calculated 


"  9T  /  p  \’ 

F8  -  (1  +  f)  ue  -  *  (1  .  f)  ue  ^  K  -fj  ,  or 

.  6  - 


f)  U 


e.eff 


> 


\ 

> 


4,  The  Cross-Sectional  Area,  A^,  at  Station  x  Along  the  Ramjet 
Engine  According  to  the  Continuity  Equation  Is 

\  ttiwi>i  nr 

Ai  ’  Ti  •  Px  V  rx'ftWlie 

For  operation  at  design  conditions  «  u^,  tWi  -  »  p^, 

and  ■  T^.  For  x  »  Dt  (diffuser  throat)  the  procedure  described  in  the 

previous  section  for  the  short  nozzle  must  be  used.  For  estimated  values 
of  and  P]Jt 


the  rji  of  the  dissociated  air  are  calculated  as  given  in  the  first  section. 
Then  the  pressure  p^.  which  is  compatible  with  the  estimated  temperature 
is  obtained  by  iteration 
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